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From a sample of 232 million T(4S') -^ BB events collected with the BABAR detector at the PEP-II B Factory in 
1999-2004, we measure the B~ -^ D^ K*~ (892) decay branching fraction using events where the K*~ is reconstructed 
in the K^t:~ mode and the D° in the K-Tr+, K-tt+tt°, and K-tt+tt+tt- channels: B{B- -^ D°K*- (892)) = 
( 5.29 ± 0.30 (stat) ± 0.34 (syst) ) x lO-". 
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From a sample of 232 million T(45') -^ BB events collected with the BABAR detector at the PEP- 
II B Factory in 1999-2004, we measure the B~ — > D _ft'*~(892) decay branching fraction using 
events where the K*~ is reconstructed in the K'^tv' mode and the D^ in the K~-k'^ , K'tt'^-k", and 
K-TY+Tv+Tv' channels: B{B' -^ D"A'*"(892)) = ( 5.29 ± 0.30 (stat) ± 0.34 (syst) ) x 10"*. 

PACS numbers: 13.25.Hw, 14.40.Nd 



The decays B^ -^ D°K^*^~ [ij are of interest because 
of their relevance to the Cabibbo-Kobayashi-Maskawa 
(CKM) model [2| of quark-flavor mixing. Interfer- 
ence effects in specific D^ final states offer a means 
of observing direct CP violation governed by the angle 
7 =arg(-KdV„;/KdV,;) i, where V is the CKM ma- 
trix. One way to access 7 is to compare the B~ -^ 
D^K*~ branching fraction to the CP-averaged branching 
fraction for a B~ to decay into a D°K*~ where the 1)° 
decays into a CP eigenstate \^. Thus a precise determi- 
nation of the B^ -^ D^K*~ branching fraction provides 
the reference for direct CP violation measurements. 

The decay B~ -^ D'^K*~ was first observed by 
CLEO 0, and later by BABAR 6]. In this paper we 
present a new measurement of the branching fraction 
B{B~ -^ D°K*~) obtained with 2.7 times more data 
than used for the previous BABAR measurement. 

This analysis uses data collected with the BABAR de- 
tector at the PEP-II e+e~ storage ring. The data cor- 
responds to an integrated luminosity of 211 fb~^ at the 
T{4:S) peak (232 million BB pairs) and 16 fb~^ at center- 
of-mass energy 40 MeV below the resonance. 

The BABAR detector is described in detail in Q. We 
give here a brief description of the components relevant 
to this analysis. Charged-particle trajectories are mea- 
sured by a five-layer double-sided silicon vertex tracker 
(SVT) and a 40-layer drift chamber (DCH) inside a 1.5 T 
solenoid. Charged-particle identification is achieved by 
combining measurements of the light detected in a ring- 
imaging Cherenkov device (DIRC) with measurements 
of the ionization energy loss (dE/dx) measured in the 
DCH and SVT. Photons are detected in a CsI(Tl) elec- 
tromagnetic calorimeter (EMC) inside the coil. We use 
GEANT4 I3] based software to simulate the detector re- 
sponse and account for the varying beam and environ- 
mental conditions. 

To reconstruct B^ -^ D^K*^ decays we select K*^ 
candidates in the K*^ -^ ^s'^^ mode and D° candi- 
dates in three decay channels: D^ — > K^tt'^ , K^tt^tt'^, 
and K~Tr~^TT~^Tr~ . Our event selection follows closely the 
one reported in [9j. K^ candidates are formed from op- 
positely charged tracks assumed to be pions with a recon- 
structed invariant mass within 13 MeV/c^ (four standard 
deviations) of the known K'^ mass, m^o [l(]j. The K^ 
candidates are fitted so that their invariant mass equals 
rrixo (mass constraint). We further require their flight 



direction and distance to be consistent with a K^ com- 
ing from the interaction point. The K^ candidate's flight 
path and momentum vectors must make an acute angle 
and the flight length in the plane transverse to the beam 
must be at least three times larger than its uncertainty. 
K*~ candidates are formed from a K*^ and a charged par- 
ticle, which are required to originate from a common ver- 
tex. We select K*~ candidates which have an invariant 
mass within 75 MeV/c^ of the known value [lOj . Finally, 
since the K*~ in B^ -^ D'^K*^ is polarized, we require 
the helicity angle Oh to satisfy |cos 9h\ > 0.35, where Oh 
is the angle in the K*'^ rest frame between the daughter 
pion and the parent B momentum. The helicity distribu- 
tion discriminates well between a B meson decay and an 
event from the e^e" '^ QQ il & {u, d, s, c}) continuum, 
since the former is distributed as cos^ 6h and the latter 
is almost flat. 

In order to reconstruct the 7r° of the D^ -^ K^tt^tt^ 
channel, we combine pairs of photons to form candidates 
with a total energy greater than 200 MeV and an in- 
variant mass between 125 and 145 MeV/c^. A mass- 
constrained fit is applied to the selected 7r° candidates. 
All D^ candidates are mass- and vertex-constrained. Par- 
ticle identification is required for the charged kaons. 
We select 13° candidates with an unconstrained invari- 
ant mass, JTiDO, differing from the world average mass, 
by less than 12 MeV/c^ for all channels except 
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K~TT~^TT^ where we require —29 < nioo — my^^ < 
+24 MeV/c^. To reduce combinatorial background in 
this channel, we further select candidates in the regions 
of the Dalitz plane enhanced by the ii:*-(892), i^*°(892) 
and ^+(770) resonances using amplitudes and phases 
measured by the CLEO experiment [ll|. In order to 
reduce the background from random two track combi- 
nations that have masses consistent with a £>" we also 
require, for the 13" -^ K^tt^ channel, |cos Od\ < 0.9, 
where Od is the angle in the Z)° rest frame between the 
daughter kaon and the parent B momentum. Finally, we 
perform a geometric fit on the B candidate which con- 
strains the 13°, the K^, and the charged pion from the 
K*~ to originate from a single vertex. 

To suppress continuum background we require 
|cos 6*^1 < 0.9, where Og is defined as the angle between 
the B candidate momentum in the T(4iS') rest frame and 
the beam axis. The distribution in cos Og is flat for qq 
events, while for B mesons it follows a sin^ 9^ distribu- 



tion. We also use global event shape variables to distin- 
guish between qq continuum events which have a two-jet 
topology in the T{AS) rest frame and BB events which 
are more spherical. We require |cos 91^\ < 0.9 where 9^ 
is the angle between the thrust axes of the B candidate 
and that of the rest of the event. We construct a lin- 
ear (Fisher) discriminant [l^l from cos 9^ and the Lq, 
L2 monomials (see below) describing the energy flow in 
the rest of the event, as in [lj|. In the center-of-mass 
frame (CM) we define Lj — Sip* \ cos 9* \^ , where i in- 
dexes the charged and neutral particles in the event once 
those from the B candidate are removed, and 9* is the 
angle of the CM-momcntum p* with the thrust axis of 
the B meson candidate. 

We identify B candidates using two nearly indepen- 
dent kinematic variables: the beam-energy-substituted 
mass TTiES = \/(s/2 -I- po • Pb^/Eq — pg and the en- 
ergy difference AE = E*^ — v^/2, where E and p are 
energy and momentum, the subscripts and B refer to 
the e"'"e~-beam-system and the B candidate in the lab 
frame, respectively; s is the square of the CM energy, 
and the asterisk labels the CM frame. 

In those events where we find more than one acceptable 
B candidate (less than 25% of selected events depending 
on the Z?° mode), we choose the one with the smallest x^ 
formed from the differences of the measured and world 
average -D° and K*~ masses scaled by the mass resolu- 
tion which includes the experimental resolution and, for 
the K*~ , its natural width. Simulations show that no 
bias is introduced by this choice and the correct candi- 
date is picked at least 80% of the time. According to 
simulation of signal events, the total reconstruction effi- 
ciencies are: 13.3%, 4.6%, and 9.0% for the D° -^ K-tt+, 
K~Tr^Tr^, and K~TT^TT^Tr~ modes, respectively. 

To study BB backgrounds we look at sideband regions 
away from the signal region in AE and m]jo. The AE 
distributions are centered around zero for signal with a 
resolution between 11 and 13 MeV for all three channels. 
We define a signal region \AE\ < 25 McV. We also define 
a AE sideband in the intervals -100 < AE < -60 MeV 
and 60< AE < 200 MeV. The lower Hmit (-100 MeV) 
is chosen to avoid selecting a region of high background 
coming from B~ -^ D*K*^. In this AE sideband we 
see no significant evidence of a background peaking near 
the B mass in niES which could leak into the signal re- 
gion. The sideband region in tti^o is defined by requir- 
ing that this quantity differs from the Z?° mass peak by 
more than four standard deviations. It provides sensi- 
tivity to doubly-peaking background sources that mimic 
signal both in AE and toes- This pollution comes from 
either charmed or charmless B meson decays that do not 
contain a true D'^. Since many of the possible contribu- 
tions to this background are not well known, we attempt 
to measure its size by including the mjjo sideband in the 
fit described below. 

An unbinned extended maximum likelihood fit to toes 



distributions in the range 5.2 < toes < 5.3 GeV/c^ is 
used to determine the event yields. For signal modes, the 
Toes distributions are described by a Gaussian function 
Q centered at the B mass with resolution (a), averaged 
over the three Z?° decay modes, of 2.7 MeV/c^. For each 
-D" decay mode k (=1, 2, 3) we determine the mean and 
sigma of the Gaussian Qk by fitting to the data. The com- 
binatorial background in the toes distribution is modeled 
with a threshold function Ak [iJl. Its shape is governed 
by one parameter ^k that is left free in the fit for each D'^ 
decay mode. We fit simultaneously tties distributions of 
nine samples: the K^ir^ , the K^tt^tt'^ and K^tt^tt^tt^ 
samples for (i) the AE signal region, (ii) the rriDO side- 
band and {Hi) the AE sideband. We fit three probability 
density functions (PDF) weighted by the unknown event 
yields. For the AE sideband, we use Ak- For the to^jo 
sideband we use N^^^p ■ Ak + N^p ■ Gk, where Qk accounts 
for the doubly-peaking B decays. For the signal region 
PDF we use N^^ ■ Ak + kN^p ■ Gk + N^ig ■ Gk, where k is 
the ratio of the mjjo signal-window to sideband widths 
and -/Vg'j is the number of B^ -^ D^ K*~ signal events. 
The AE sideband sample helps define the shape of the 
background function Ak- We assume that the B decays 
found in the mjjo sideband have the same final states 
as the signal so we use the same Gaussian shape for the 
doubly-peaking B background. 



TABLE I: Results from the fit and quantities used to derive 
the B~ — > D^K*~ branching fraction. For each channel we 
give the event yield resulting from the fit, the efficiency, and 
the branching fraction measurement, in units of 10~*, derived 
using Eq. Q The uncertainties are statistical only. 



K 



K- 



K- 



TT 



Yield 
Efficiency 
B{B- -^ D°K'- 



144 ± 13 185 ± 19 

13.30% 4.60% 

5.15±0.47 5.65±0.57 



195 ± 18 

8.99% 
5.24±0.49 



The fit results are shown graphically in Fig. Q] and nu- 
merically in Tabled For each channel k, a measurement 
Bk of the branching fraction B{B^ -^ D'^K*~) is derived 
as follows: 



Bk^ 



NiD" -^Xk)-f 



Ngi-ek-BK- -BiD^^Xk)' 



(1) 



where N{D^ -^ Xk) is the event yield from the fit, / the 
fraction of K*~'s in the sample (discussed below), Nb± 
is the number of charged B mesons in the data sample, 
£fe is the efficiency to reconstruct B~ -^ D^K*~ when 
D° -^ Xk, Bk- = B{K*- -^ K^tt) ■ B{K° -^ tt+tt") 
and B{D^ — > Xk) are the branching fractions of the K*^ 
and the D^ . We have assumed equal production of pairs 
of neutral and charged B mesons in T{4,S) decay. 

Systematic effects arise from the difference between the 
actual detector response for the data and the simulation 
model for the Monte Carlo. Here the main effects stem 
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When we remove the requirement on the K*~ hehcity an- 
gle, we see that the K*~ hehcity distribution (Fig.jSJ of 
the selected events manifests a forward-backward asym- 
metry that indicates an interference with a K'^tt back- 
ground [3, llal- We model the K^'!t~ system with a P- 
wave and an S-wave component. The P-wave mass de- 
pendence is described by a relativistic Breit-Wigner while 
the S-wave piece is assumed to be a complex constant. 
This model is fitted to the data and shown in Fig. El 
along with an estimate of the combinatorial background. 
Neglecting higher resonances, the number of K'^tt~ peak- 
ing background events is (4 ± 1)% of the total measured 
number of signal events. We do not quote a systematic 
error on the contributions of the neglected partial waves 
{non-K* P-wave and higher order waves) since their ex- 
pected rates in the Kt: mass window are far below that of 
the S'-wave ^^. In Fig. 13 we see that a relativistic Breit- 
Wigner gives a fair description of the resonance structure 
in the K'^tt~ mass spectrum (x^=26.8 for 20 degrees of 
freedom) . 



FIG. 1: Distributions of m-ES in the signal region for B -^ 
D°K*- decays where D° -^ K--R+ (top), A'"7r+7r° (mid- 
dle), and K~-K^-K~-K^ (bottom). The dashed curve indicates 
the contribution from the combinatorial background and the 
peaking B-background which is estimated from a simultane- 
ous fit to the D" sideband (not shown). 



from the modeling of the tracking efficiency (1.2-1.3% 
per track), the K^ reconstruction efficiency (2% per K^)^ 
the tt" reconstruction efficiency for the K~-K^'n^ channel 
(3%) and the efficiency and misidentification probabil- 
ities from the particle identification (2% per kaon). A 
study of a high-statistics B~ -^ D°Tr~ control sample 
shows excellent agreement between the data and Monte 
Carlo sample except for the distributions of AE and the 
continuum-suppression Fisher discriminant. For these 
variables, differences of up to (2.5 ± 1.1)% are measured 
between the data and Monte Carlo. Suitable corrections 
to the efficiencies are therefore applied and systematic er- 
rors assigned. The K*~ hehcity angle distributions differ 
significantly between data and simulation because of the 
non-resonant background under the K*~ peak. We de- 
scribe below how we subtract this background. For the 
pure K*~ events, we estimate that the residual discrep- 
ancy between data and simulation in the helicity to be 
less than 1.6%. We determine using simulations that the 
TTiEs signal PDFs deviate from the single Gaussian shape 
by less than 0.1%. Substantial systematic uncertainties 
come from the measured D'^ branching fractions [10| and 
the number of B pairs in the sample. 

The observed number of signal events must be cor- 
rected for the non-resonant K'^tt pairs under the K*~ . 



TABLE II: Systematic uncertainties. Xk refers 
decay modes given in the columns. S^.- is the 



fraction of the observed K* 
chain. 



Kln- 



m 



to the D° 
branching 
^7r~ decay 



Source 


K-n+ 


_ft:-7r+7r" K-7T+7T-n+ 


Tracking efficiency 


3.8% 


3.8% 


6.3% 


tt" efficiency 


- 


3.1% 


- 


Particle Identification 


2.0% 


2.0% 


2.0% 


K^ efficiency 


1.6% 


1.9% 


1.8% 


coseniK'-) 


1.6% 


1.6% 


1.6% 


Fisher 


1.1% 


1.1% 


1.1% 


AB 


1.9% 


1.8% 


2.0% 


rriES PDF shape 


0.1% 


0.1% 


0.1% 


Number of B^ 


1.1% 


1.1% 


1.1% 


Simulation statistics 


0.9% 


1.4% 


1.0% 


B^,. [10] 


0.2% 


0.2% 


0.2% 


B(d0 ^ Xfe) [lOj 


2.4% 


6.2% 


4.2% 


K^ TT~ S-wave subtraction 


1.1% 


1.1% 


1.1% 


Total systematic error 


6.1% 


9.0% 


8.7% 



All sources of systematic uncertainties are listed for 
each mode in Table E] With the exception of AE and 
simulation statistics the systematic error sources listed 
in Tableware correlated among the different Z?° modes. 
We use the procedure discussed in llfl| to form a weighted 
average of the three D*^ decay modes and determine: 



B{B- -^ D°K*-) = ( 5.29 ± 0.30 ± 0.34 



X 10" 



The first error is statistical and the second is systematic. 
We have compared the results from this analysis using 
the same data set as in our previously published analy- 
sis [^. The two analyses use different selection criteria 
and therefore find different numbers of events. The re- 
sults from the two analyses are consistent to within a 
half of a (statistical) standard deviation. We have also 
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FIG. 2: Acceptance corrected distribution of cos 6h{K*~). 
The solid line is a fit to a model which includes P-wave and 
S-wave interference. The dotted line shows the combinatorial 
background as estimated from the data. 
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FIG. 3: Invariant mass of KgTT~ combinations with all other 
analysis cuts applied. The solid curve is a Breit-Wigner line 
shape including detector resolution. The dotted line shows 
the combinatorial background. 



calculated the branching fraction for the two data sets 
obtained since the previous analysis. The measurement 
in each set is consistent with, although lower than the 
value obtained in |y|. This result supersedes our previ- 
ously published result. 

In summary, we have measured the branching fraction 
of the decay B- -^ D"K*- in the D°K°Tr- final state 
and observed the interference of the K*~ with a small 
non-resonant K'^tt background. 
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